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Dark energy constraints have forced viable alternatives that differ substantially from a cosmo- 
logical constant A to have an equation of state w that evolves across the phantom divide set by 
A. Naively, crossing this divide makes the dark energy gravitationally unstable, a problem that is 
typically finessed by unphysically ignoring the perturbations. While this procedure does not affect 
constraints near the favored cosmological constant model it can artificially enhance the confidence 
with which alternate models are rejected. Similar to the general problem of stability for w < 0, the 
solution lies in the internal degrees of freedom in the dark energy sector. We explicitly show how 
to construct a two scalar field model that crosses the phantom divide and mimics the single field 
behavior on either side to substantially better than 1% in all observables. It is representative of 
models where the internal degrees of freedom keep the dark energy smooth out to the horizon scale 
independently of the equation of state. 



I. INTRODUCTION 

A self-consistent model for the dark energy requires not 
only a parameterization of the evolution of its equation 
of state in the background w — p/p but also a physical 
model for its spatial fluctuations to guarantee gravita- 
tional stability. Cosmological constraints on a constant 
equation of state have continued to close in upon w = — 1 
(e.g. 0). Since models with w < — 1 are called phantom 
dark energy models (e.g. Q), we call this the "phantom 
divide" . Viable alternate models for a strongly evolving 
w therefore must cross the divide at intermediate redshift 
so that the effects on either side cancel. Simple general- 
izations of a single scalar field model that cross the divide 
in fact cause severe gravitational instabilities in the dark 
energy sector (e.g. |3|). 

The usual approach in the literature for dealing with 
such cases is to artificially turn off the dark energy per- 
turbations. Doing so violates energy-momentum conser- 
vation whenever w ^ — 1. The justification for dropping 
these perturbations is that observations already place w 
close to —1 and so the induced error is presumably small 
in some physical realization of a crossing model. While 
true for the currently allowed la deviations of Aw ~ 0.1, 
the confidence level at which larger deviations can be re- 
jected can be affected. Furthermore, a strong time evolu- 
tion allows w to differ substantially from — 1 during some 
epochs and still be consistent with the distance data, ft 
is therefore important to show explicitly that models ex- 
ist where the dark energy remains smooth as it crosses 
the phantom divide as implied by the usual procedure. 

The need for a self-consistent treatment of the evolu- 
tion of the dark energy is most apparent for the cosmic 
microwave background. Here the ISW effect is sensitive 
to the decay of the gravitational potential and, for ex- 
ample, the difference between the smooth and clustered 
regimes of the dark energy for a constant w w —2/3 is 
roughly 50%. While the impact for a canonical scalar 
field which is smooth out to the horizon scale is lower, 
it is well known that it changes CMB predictions sig- 



nificantly for larger w and hence the confidence level of 
constraints on highly deviant w (e.g. 

In this Brief Report, we explicitly construct a two 
scalar field model of the dark energy that is gravitation- 
ally stable across the phantom divide and matches the 
canonical single scalar field predictions on either side to 
much better than 1% in all observables. Taken literally, 
such a model of course compounds the coincidence prob- 
lem of dark energy but as we shall show it is a proxy 
for a potentially wider class of dark energy models whose 
internal degrees of freedom keep it smooth out to the 
horizon scale. 



II. INSTABILITY 

It is well known that dark energy models beyond a 
cosmological constant require internal degrees of free- 
dom, or the presence of non-adiabatic stress perturba- 
tions, to remain gravitationally stable. This necessity 
arises from momentum conservation. Consider the di- 
mensionless momentum density pui = T t of the dark en- 
ergy stress tensor TJ\,. The scalar component in Fourier 
space (e.g. [f| and 0| for a pedagogical treatment in the 
same notation) 

pii = — (3w — l)pu + kSp + (1 + w)pkA, (1) 
a 

where 2 A = 5g 00 /g Q is the time-time perturbation to 
the metric in an arbitrary gauge. Given dark energy 
fluctuations that are internally adiabatic p(p) (not to be 
confused with adiabatic across all energy density compo- 
nents) 

fp\ ( l d\n{l + w) \ 

where w = p/p in the background. Perturbations go un- 
stable whenever the pressure response to a density fluctu- 
ation is negative or singular. The former occurs for adi- 
abatic pressure fluctuations if w < and 1 + w is slowly 
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varying. A viable dark energy candidate in this regime 
must contain internal degrees of freedom that supply a 
non adiabatic pressure. For a single scalar field with a 
canonical or non-canonical kinetic term, this is achieved 
through separate kinetic and potential contributions to 
the energy density and pressure. 

In this case, one defines a sound speed c 2 which re- 
lates the energy density and pressure fluctuation of the 
kinetic term |8j or equivalently of the zero momentum or 
constant field gauge and obtains 
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If the sound speed c e < 1 it determines the scale under 
which the dark energy is effectively smooth through the 
sound horizon J c e da / (a 2 H) . For the canonical kinetic 
term c e = 1. More specifically, well above this scale, 
stress gradients are negligible and the gravitational po- 
tential A = ^ in the Newtonian gauge evolves as 
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Well below this scale the dark energy is smooth and VP oc 
G with 
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+-[l-i£7(a)]n DE (a)G = 0, (4) 

where f2DE(a) = 8irGp/3H(a) 2 and a = 1 is assumed 
where no dependence is given. Note that in both limits, 
the effect of the dark energy on the potential is solely 
a function of its background energy density. The true 
degree of freedom in a dark energy model is where this 
transition occurs. Any physical solution to the instabil- 
ity problem that matches a desired w(a) and transition 
scale will be a fairly robust representation of the class of 
models. We will use this fact below to replace the usual 
single scalar field ansatz with two scalar fields. 

A single scalar field does not generally solve the prob- 
lem that p/p becomes singular as w evolves across 
the phantom divide w = — 1 with finite slope d(l + 
w)/d\na ^ since p/p still appears in Eqn. @. Sta- 
bility can obviously be achieved by an alternate ansatz 
for the internal degrees of freedom. The simplest solution 
that preserves the behavior of the single scalar field tran- 
sition scale away from the crossing point is to introduce 
multiple scalar fields. 



III. TWO FIELD MODEL 

For definiteness, the target form of w(a) that we wish 
to model with two scalar fields is 



w(a) = w + (1 - a)w a , 



(5) 



where w a is a constant. For simplicity, we will take the 
two fields, denoted "+" and "— " to individually have 
constant equations of state 



w± = w ± e r , 



(6) 



Note that this restriction to strictly constant equations 
of state is not essential to the construction. Though a 
strictly constant w for a rolling scalar field is unrealistic, 
scalar field potentials do exist where the resultant equa- 
tion of state differs significantly from ±1 and roughly con- 
stant during the redshifts of interest [lj} ■ I n anv case the 
point of this explicit construction is to provide a cross- 
ing model that is simple to implement in existing cos- 
mological codes and does not violate energy-momentum 
conservation. It is not intended to be a well-motivated 
model. 

The equations of state define the relative energy den- 
sity contributions as a function of redshift or scale factor 
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Here S n — S(a n ) defines the ratio at a normalization 
epoch a n . This epoch should not be chosen as a n = 1 
since variations in the equation of state locally leave no 
net effect. Rather a n = 3/4 is roughly the pivot point 
where variations in the equation of state make the max- 
imal effect on the high and low redshift observables. We 
will adopt this value as the matching point between the 
two field model and the target w(a). 

With equal sound speeds in the two components, the 
pressure fluctuation becomes 



Sp = c 2 5p + 3-- \c 2 pu 
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w + p + u + — w-p-uJ\ , (9) 



Here p = p+ + P-,P = P++P-, pu = p+u+ + p-U- none 
of which contain singularities at the crossing. 

The model dark energy equation of state becomes 



ad (a) = w + S(a)e n 



(10) 



and at the normalization point w n — w + S n e n . Likewise 
at the normalization point the derivative of w is 
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Thus w a > if S 2 < 1. A negative derivative is possible 
if Si > 1, i.e. if the p_ < 0. 

The two field model has 3 parameters w, e„ and S n . 
The target model has 2 parameters wq and w a , leaving 
one adjustable parameter to improve the performance of 
the parameterization. Choosing this parameter to be e„ 
defines the other two in terms of wq, w a and e„ 
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wq + (1 - a n )w a - S n e n , 



(12) 
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and S„ 
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We now choose e„ to satisfy several criteria. Firstly, as 
w a — > the two field model should reduce to a single field 
so as to exactly match the standard dynamics. Secondly, 
the two field model should match a sufficiently wide range 
of -Wm&x <w a < w max . We will take ui max = 1. Finally, 
it should minimize higher order derivatives. 

To satisfy the first condition, e n should scale as a power 
of w a at least near w a — 0. Combined with the second 
criteria. 
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(Vj a > 0) (14) 



since w a is maximized at S n = 0. 

The final condition is that higher order derivatives such 

as 
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should be minimized in so far as possible. At large w a 
this is ensured by keeping S n e n small. A good choice to 
minimize these higher order derivatives is to take p = 
1/10. For w a < 0, S n > 1 and so one must restrict the 
values of e„. We take 



max(e„,0.1) , (w a < 0) . 



(16) 



Since 1 — 5^ is unbounded from below, this restriction still 
allows large negative w a in Eqn. (|lll) while preserving the 
single field correspondence for w a — * from below. 

To complete the modelling we match the physical en- 
ergy density (i.e. not relative to critical) at the normal- 
ization epoch a n = 3/4. To scale the dark energy den- 
sity of the target dark energy model Eqn. JSJ we write 
Pde(o„) = ^DE/0 C rit5DE where 



3DE = a - 3{1+W0+w ^e' s{1 - a " )% 



and analogously introduce 



g± = a" 3 ^^) 



(17) 
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Thus given a target model with an energy density relative 
to critical of Ode and Hubble constant parameterized 
by Ho = 100/ikm s _1 Mpc -1 , the 2 field model has an 
effective Hubble constant h c 
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+ 1-«de (19) 



and density relative to critical defined by that Hubble 
constant of 



fl± — f^DEffDE" 



1±6„ 
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(20) 
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FIG. 1: Target versus 2 field model behavior of w and the 
dark energy observables relative to a fiducial A model (dotted 
line, see text). Even in this nearly worst case scenario where 
w crosses the phantom divide near the pivot point a ~ 3/4 
or z ~ 0.3, the model reproduces the observables to 0.1% 
accuracy for w a > and only a factor of a few larger for 
w a < 0. 



Although this construction is completely phenomeno- 
logical and hence physically contrived, the general point 
is that once a close matching of w(a) to some target has 
been achieved with multiple scalar fields, the scalar field 
dynamics will make the predictions in both the smooth 
and clustered regimes robust to reparameterization. 



IV. DISCUSSION 

In Fig. n we show examples of the performance of 
the model on the dark energy observables of comov- 
ing distance D(z), Hubble parameter H(z) and linear 
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FIG. 2: Model predictions for the CMB anisotropy power 
spectrum relative to a fiducial ACDM model. For models 
that do not cross the phantom divide, e.g. wo = —I, w = 0.5, 
the 2 field model matches the single field predictions to 0.1%. 
For crossing models that mimic a cosmological constant by 
having w n ~ — 1, e.g. wo = —1.15, w = 0.5, the deviations 
are essentially indistinguishable. 

growth rate G(z) in the smooth dark energy regime (see 
Eqn. for strong variations in the equation of state 
w a = ±0.5. To better show the relevant performance 
we choose wq such that w n « —1 and plot the observ- 
ables relative to a nearly degenerate fiducial A model of 
^de = 1 — &m = 0.74 and h = 0.735 which is a good 
fit to the current CMB data. In all cases we adjust the 
firjE of the target model to match the comoving angular 
diameter distance to recombination and hence the CMB 
peak results. Note that for w a > 0, despite a mismatch 
in w at very high redshift, all observables as a function 
of redshift are modeled to the 0.1% level. The perfor- 
mance for w a < is somewhat worse but never exceeds 
a few times this level. Thus the two field model will 
remain an adequate parameterization until the statisti- 
cal and systematic errors in the dark energy observable 
measurements reaches the sub 1% level. Mismatches at 
this point merely reflect the unavoidable fact that higher 



order derivatives in w(a) produce observable effects and 
our target constant w a model is itself inadequate. Even 
in this regime the 2 field model can be useful since the 
parameter e„ can be used to marginalize or probe the 
second derivative. 

The two field model allows one to also calculate the 
CMB anisotropy. For a case with no crossing of the di- 
vide, e.g. w a — 0.5, wo — — 1, the single field model with 
the exact w(a) matches the 2 field model with the ap- 
proximate w(a) to ~ 0.1%, in particular at the low mul- 
tipoles of the ISW effect (see Fig- EJ) - In addition to the 
adequate matching of w(a), this indistinguishability is a 
consequence of setting all of the sound speeds to c e = 1. 
Aside from w(a) the remaining degrees of freedom in a 
dark energy model involve the transition scale between 
the smooth and clustered regimes. We also show the pre- 
dictions for a model that does cross the phantom divide 
w a = 0.5, wo = —1.15. As this model has an equation of 
state at the normalization point of w n = —1.025, it pre- 
dicts even smaller deviations from the fiducial w = — 1 
model with no instabilities in the evolution of the gravi- 
tational potential. 

The model constructed here permits a self-consistent 
likelihood analysis of dark energy observables involving 
both models that cross the phantom divide and those 
that differ strongly from A in the pivot w n . Although pre- 
dictions for the former class models differ little from the 
currently favored smooth cosmological constant case, the 
confidence level at which the latter can be excluded can 
be affected by a consistent model of dark energy cluster- 
ing. A detailed study of the effect on current cosmological 
constraints is beyond the scope of this Brief Report. 

If future observations require an evolutionary crossing 
of the phantom divide, it will be a good indication that 
the dark energy contains hidden internal degrees of free- 
dom in its physical structure. 
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Note added in proof: After completion of this work, 
we became aware of related work by Guo et al 
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